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Abstract 

In this comment]] I show that the experimental data on quantum diffusion 
of 3 He impurities in solid 4 He can be explained using the adopted quasipar- 
ticle theory. The contention by E.G. Kisvarsanyi and N.S. Sullivan (KS) in 
Phys.Rev.B 48 16557 (1993) as well as in their Reply (ibid. 55 3989 (1997)) 
to the Grigor'ev's Comment (Phys.Rev.B 55 3987 (1997)) that "Pushkarov's 
theory of phonon scattering fails to fit the data by very large factors" is 
groundless and may result from their bad arithmetical error. This means 
that the phonon-impurity scattering mechanism of diffusion is consistent with 
experiment and its neglecting by KS makes their results questionable. 
The temperature dependence of the impuriton diffusion coefficient in solid helium can 
be determined by two types of thermal excitations - phonons and vacancies. Both of them 
have been considered in a number of works (see e.g. |l]Q]). For low concentrations and tem- 
peratures T < 1 K the theory based on the impuriton-phonon scattering has been confirmed 
by the experiment (cf e.g. |1]|||3J] anc ^ the references therein) and used for determination of 
quasiparticle characteristics. The quantitative agreement was first obtained in Ref. || (cf 
also ||[7|]), an d since that time the theory has not been subjected to principle changes. 

Kisvarsanyi and Sullivan (KS) argue in their work || and in their Reply || to the 
Grigor'ev's Comment |TD| that they have proposed ,,a new theoretical treatment of the 



temperature dependence of the diffusion of isotopic impurities in solid 4 He" as well as that my 
theory of phonon scattering ,, fails to fit the data by very large factors". The latter assertion 
is very important because it concerns the generally accepted selfconsistent approach to the 
diffusion in quantum crystals (see e.g. [Jl||12],|]-0,|^-|l^,|l|,^J5]). It has been used as the only 
argument to neglect fully the phonon-impurity scattering mechanism, and turn back to the 
vacancy controlled impurity diffusion. 

KS have evaluated in || the factor A in the temperature dependence of the diffusion 
coefficient D = AT~ 9 . They have found A "to be in the range A ~ 10 4 — 10 5 " while 
my theoretical prediction p|-f^pl2|-|14|l is A ~ 10 7 cm 2 s' 1 K~ 9 and the experimental values 



are A « 10 7 (Ref. 0) or A « 2.4 x 10 7 (Ref. §]). They concluded that "the phonon 
scattering is too weak by a factor of at least 100 to explain the observed diffusion" 



1 This Comment was sent to the Editor of Phys.Rev.B on 03.12.1998, accepted for publication in 
August 1999, but not published till now. 
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p. 16579). KS have used in their calculation eqs. (15) and (16) given below as (1) and (2) 
respectively: 

D = \za 2 J 2 T « 4.3 x 1(T 15 J 2 r, (1) 



D * uj d 1 ^ 0.79 x 10- 17 e^T- y (2) 



where 9^ is the Debye temperature, ujd = kQ D /h, a = 3.27 x 10 8 cm is the interatomic 
distance, J is the tunneling frequency and s — 1/3. Hence, 

A = 3.4 x 10~ 32 J 2 6^ (3) 

Although these expressions do not exactly coincide with those used in my works they should 
give a correct order of magnitude for the diffusion coefficient and, in particular, the order of 
A. KS have obtained for Q D = 30K and J = 2.5 MHz the value A = 6.0 x 10~ 5 crrfs^T 9 
while the right value (obtained by substituting the same numbers into (3)) is A = 1.4 x 10 _r . 
The issue is not that the above expressions are the most fundamental or give the best value 
of A, but rather that they yield A values consistent with experiment. The value obtained 
by KS and leading to the rejection of the phonon scattering mechanism (cf M, p. 16579) 
may come, therefore, from an arithmetical error. 

In fact, 0£> = 30 K is the upper limit and corresponds to a molar volume V m = 19.8 
cm 3 while the experimental values cited above are for V m = 21 cm 3 with = 26 K 
Hj| . This yields A = 0.44 x 10 T . An ambiguity appears with the notation J = 2.5 



MHz (p. 16579) and J 34 /27r = 2.3 x 10 5 s" 1 = 0.23 MHz (p. 16580) used for one and 
the same quantity. I suppose that J is given in rad/s. Otherwise, it should have been 
written in the form J/2tt = 2.5 MHz (as for J 34 ) and would be more than 10 times larger 
than J 34 leading to J = 2vr2.5 x W 6 h/k B = 1.2 x 10" 4 K and to an energy band width 
A = 2zJ »3x 10~ 3 K in a drastic disagreement with all known experiments. On the other 
hand the value J34 ~ 1.10 x 10 -5 K is typical of the exchange integral of 3 He atoms in solid 
4 He and its substitution into (3) gives the correct order of magnitude for A. Finally, even 
if the unrealistic value J = 2n2.5 ~ 15.7 MHz were used in (3) the result for A differs from 
that of KS by more than an order of magnitude. 

It is clearly seen, therefore, that expressions (1) and (2) (not necessarily the best ones) 
are in agreement with the experimental data. This can be easily verified by substituting 
J and ©£. Therefore, the argument that ,,the diffusion constant calculated for this theory 
fails to fit the experimental data by a factor of 100" || fails. There are no experimental 
data to require any drastic change of the eqs. (1) and (2). 

It is not accurate to cite my paper || in [|J (Ref.18) as if it were in support to the values 
of A ,,in the range 10 -4 — 10~ 5 ". The corresponding value is A = 2 x 10 -7 as follows from 
eqs.(l) and (10) in ||. It was the first correct evaluation of the impuriton band width. It is 
worth noting in addition that the idea of vacancy controlled mechanism could not be called 
new because it was considered about 20 years ago (cf e.g. 0) and was found not satisfactory 
(cf e.g. JTJ) to explain the temperature dependence of the diffusion coefficient. 

As a consequence, the good fit of the vacancy mechanism reported by KS after neglecting 
the phonon-impuriton scattering gets at least doubtful. There are also other circumstances 
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which call the impurity diffusion description by KS || in question, some of them being listed 
in |To| (see also ||16|| ). They will be considered elsewhere. 

I do not concern the problem of impurity diffusion in solid hydrogen as not relevant to 
the concrete discussion. If the theory of KS does not work well for helium, it obviously 
cannot be a ,, universal theory" for both quantum solids as argued in ||. 

In addition, the work of KS || suffers from a number of inaccuracies. If the dispersion law 
is defined by their eq.(3), then the energy band width in a simple lattice is A = 2zJ, not zJ, 
the circles and the squares in Fig. 1 have to change places, the works of Landesman (their 
[13]) are published in Ann. de Phys. (French), not in Ann. Phys. (N.Y.), the number 1/3 in 
eq.(5) for the dispersion law in a bcc lattice has to be replaced by 1/2, the term cos (^£^a ) 

in (6) has to be replaced by cos (^/|&; 2 ao). The correct formulae can be found e.g. in my 
works []1^,|6|) |15| ) |I4[ . The problem is not only in misprints - KS present as their own results 
known for 25 years. The velocity squared < v 2 >= l&a^J 2 (eq.(7) in [||) was calculated 
first by Sacco and Widom |18[ in 1976, not in their work |19| of 1992 as cited (it is another 
question whether it is the right quantity for the problem under consideration). The way of 
calculation of vacancy-impurity cross-section <7o = 1.40a 2 is a secret. Following the references 
one sees that one of the authors has made a private communication to herself (Ref.13 in the 
work of Kisvarsanyi, Runge and Sullivan |2(J reads "K.Runge, private communication", and 
Ref. 14. is an unpublished M.S. Thesis of Kisvarsanyi). 



3 



REFERENCES 



[1] A. R. Allen, M. G. Richards and J. R. Schratter, J. Low Temp. Phys. 47 289 (1982) 
[2] D. R Locke, J.Low Temp.Phys. 32, 159 (1978) 

[3] V. A. Mikheev, B. N. Essel'son, V. N. Grigor'ev and N. P. Mikhin, Fiz.Nizkih Temp. 3 
385 (1977) 

[4] B. N. Essel'son et al, in Proc. 20-th Soviet Conf. on Low.Temp.Phys. NT-20, v.2 (1979) 

p.236; J.de Phys.Colloque C6, suppl.8, 39 119 (1978) 
[5] D. I. Pushkarov, Pisma Zh.Eksp.Teor.Fiz. 19, 751 (1974) [JETP Lett. 19(12) 386-387 

(1974)] 

[6] D. I. Pushkarov, Fiz.Nizkih Temp. 1 581 (1975); ibid 1 585 (1975) 
[7] D. I. Pushkarov, Zh.Eksp.Teor.Fiz. 68 1471 (1975); [Sov.Phys. JETP 41 (4) 735-737 
(1976)] 

[8] E. G. Kisvarsanyi and N. S. Sullivan, Phys.Rev.B 48 16557 (1993) 
[9] E. G. Kisvarsanyi and N. S. Sullivan, Phys.Rev.B 55 3989 (1997) 
[10] V. N. Grigor'ev, Phys.Rev.B 55 3987 (1997) 

[11] D. I. Pushkarov, Zh.Eksp.Teor.Fiz. 59, 1755 (1970) [Sov.Phys. JETP 32(5) 954 (1971)] 

[12] D. I. Pushkarov - "Quantum Theory of Defects in Crystals at Low Temperatures", PhD 
Thesis, Moscow University, Moscow, 1972 (in Russian). 

[13] D. I. Pushkarov, Quasiparticle Theory of Defects in Solids, World Scientific, Singapore- 
New Jersey-London-Hong Kong, 1991. 

[14] D. I. Pushkarov, Defektony v Kristallakh . . . (Defectons in Crystals. Quasiparticle Ap- 
proach to Defects in Solids), Moscow, Nauka 1993. (in Russian.); [first edition, JINR 
PI 7-87- 177 (1987), Dubna, Russia]. 

[15] D. I. Pushkarov, Quasiparticle Approach in Quantum Theory of Solids, DSc Thesis, 
Dubna, USSR, 1986. (in Russian.); Preprint JINR 17-86-223 (1986), Dubna, Russia 

[16] V. N. Grigor'ev, Fiz.Nizkih Temp. 23 5 (1997) 

[17] D. O. Edwards and R. C. Pandorf, Phys.Rev. 140, A816 (1965) 

[18] J. E. Sacco, A. Widom, J.Low Temp.Phys. 24 241 (1976) 

[19] M.Rall, D.Zhou, E.K.Kisvarsanyi and N.S.Sullivan, Phys.Rev.B45, 2800 (1992) 

[20] E. G. Kisvarsanyi, K. Runge and N. S. Sullivan, Phys.Lett.A 155 337 (1991) 



4 



